eluding the basic principles that a starting 
configuration is randomly selected from 
within the parameter space, the algorithm 
tests other configurations with the goal of 
finding the globally optimal solution, and 
the region from which new configura- 
tions can be selected shrinks as the search 
continues. The key difference between 
these algorithms is that in the SA algo- 
rithm, a single path, or trajectory, is taken 
in parameter space, from the starting 
point to the globally optimal solution, 
while in the RBSA algorithm, many trajec- 
tories are taken; by exploring multiple re- 
gions of the parameter space simultane- 


ously, the algorithm has been shown to 
converge on the globally optimal solution 
about an order of magnitude faster than 
when using conventional algorithms. 

Novel features of the RBSA algorithm 
include: 

1. More efficient searching of the pa- 
rameter space due to the branching 
structure, in which multiple random 
configurations are generated and 
multiple promising regions of the pa- 
rameter space are explored; 

2. The implementation of a trust region 
for each parameter in the parameter 
space, which provides a natural way of 


enforcing upper- and lower-bound 
constraints on the parameters; and 

3. The optional use of a constrained gra- 
dient-search optimization, performed 
on the continuous variables around 
each branch’s configuration in param- 
eter space to improve search efficiency 
by allowing for fast fine-tuning of the 
continuous variables within the trust 
region at that configuration point. 

This work was done by Matthew Bolcar, J. 
Scott Smith, and David Aronstein of Goddard 
Space Flight Center. Further information is con- 
tained in a TSP (see page 1). GSC-15908-1 


@ Method for Pre-Conditioning a Measured Surface Height Map 
for Model Validation 

This method can be implemented in most optical modeling and simulation software. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


This software allows one to up-sample 
or down-sample a measured surface 
map for model validation, not only with- 
out introducing any re-sampling errors, 
but also eliminating the existing meas- 
urement noise and measurement errors. 
Because the re-sampling of a surface 
map is accomplished based on the ana- 
lytical expressions of Zernike-polynomi- 
als and a power spectral density model, 
such re-sampling does not introduce any 
aliasing and interpolation errors as is 
done by the conventional interpolation 
and FFT-based (fast-Fourier-transform- 
based) spatial-filtering method. Also, 
this new method automatically elimi- 
nates the measurement noise and other 
measurement errors such as artificial 
discontinuity. 

The developmental cycle of an optical 
system, such as a space telescope, in- 
cludes, but is not limited to, the follow- 
ing two steps: (1) deriving requirements 
or specs on the optical quality of individ- 
ual optics before they are fabricated 
through optical modeling and simula- 
tions, and (2) validating the optical 
model using the measured surface 


height maps after all optics are fabri- 
cated. There are a number of computa- 
tional issues related to model validation, 
one of which is the “pre-conditioning” 
or pre-processing of the measured sur- 
face maps before using them in a model 
validation software tool. 

This software addresses the following 
issues: (1) up- or down-sampling a meas- 
ured surface map to match it with the 
gridded data format of a model valida- 
tion tool, and (2) eliminating the sur- 
face measurement noise or measure- 
ment errors such that the resulted 
surface height map is continuous or 
smoothly-varying. So far, the preferred 
method used for re-sampling a surface 
map is two-dimensional interpolation. 
The main problem of this method is that 
the same pixel can take different values 
when the method of interpolation is 
changed among the different methods 
such as the “nearest,” “linear,” “cubic,” 
and “spline” fitting in Matlab. The con- 
ventional, FFT-based spatial filtering 
method used to eliminate the surface 
measurement noise or measurement er- 
rors can also suffer from aliasing effects. 


During re-sampling of a surface 
map, this software preserves the low 
spatial-frequency characteristic of a 
given surface map through the use of 
Zernike-polynomial fit coefficients, 
and maintains mid- and high-spatial- 
frequency characteristics of the given 
surface map by the use of a PSD model 
derived from the two-dimensional PSD 
data of the mid- and high-spatial-fre- 
quency components of the original sur- 
face map. Because this new method 
creates the new surface map in the de- 
sired sampling format from analytical 
expressions only, it does not encounter 
any aliasing effects and does not cause 
any discontinuity in the resultant sur- 
face map. 

This work was done by Erkin Sidick of 
Caltech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1). 

The software used in this innovation is 
available for commercial licensing. Please con- 
tact Daniel Broderick of the California Insti- 
tute of Technology at danielb@caltech.edu. 
Refer to NPO-47593. 
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